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Abstract 

The mass estimate of the d*{IJ^ = 03^) dibaryon is improved by a dynam- 
ical calculation in the quark delocalization, color screening model. The partial 
decay width of d* into an NN D-wave state is also obtained. The mass ob- 
tained is slightly larger than that obtained in adiabatic calculations, due to the 
anharmonicity of the effective potential between two A's. The value of the width 
obtained due to tensor one-gluon-exchange is about 5 MeV, comparable in mag- 
nitude to earlier results found using pion exchange. 
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1 Introduction 



Quantum Chromodynamics (QCD) is widely accepted as the fundamental theory of 
the strong interaction. However, the low energy behavior of QCD remains a challenge 
to understand. Lattice QCD can calculate the structure of simple quark-gluon systems 
but it remains extremely difficult to calculate complicated systems such as two-hadron 
interactions and the structure of multiquark systems. Quark models were inspirational 
in the development of QCD and may play a similar role in understanding the low 
energy behavior of QCD. 

The nonrelativistic constituent quark model has been quite successful in understand- 
ing hadron spectroscopy This model has been extended by many groups to the 
description of baryon-baryon (B-B) interactions. There is a broad consensus that the 
repulsive core of the N-N interaction can be attributed to the quark internal structure 
of the nucleon. In this respect, the repulsive core is quite similar to that of molecular 
forces originating from the electronic internal structure of atoms. However, the inter- 
mediate range attraction of the N-N interaction remains absent in many quark model 
approaches and a scalar meson exchange has frequently been invoked to remedy this. 

Spontaneous chiral symmetry breaking supports the view that Goldstone boson ex- 
changes occur even between constituent quarks. Based on this, recently a phenomeno- 
logical quark model has been proposed by Glozman et al. [0. The model gives a rather 
good description of the baryon spectrum and it has also been applied to NN interac- 
tions and dibaryons 0].But a question remains as to what extent this picture should 
be used: Should it be used for long range pion exchange only, or extended to interme- 
diate two pion or a exchange or even to include short range vector and other heavy 
meson exchange? At what point should gluon exchange between quarks be replaced? 
Of course, a great deal of further study is needed to settle such questions definitively. 

One point we would like to emphasize here is that such an approach would qualitatively 
differentiate between nuclear and molecular forces, which have been known for more 
than half a century to be similar in many features, the energy and length scale differ- 
ences notwithstanding 0. Recently, Anderson has pointed out that it is questionable 
to attribute the N-N interaction to meson exchange P| except for the longest range 
portion of pion exchange. NN partial wave analysis can only justify long range pion 
exchange but not other, heavier meson exchanges 0. 

In the quark model approaches mentioned above, a two-body confinement potential is 
usually used. We recognize that this is a highly simplified model assumption. Quark 
confinement is a nonperturbative property of low energy QCD. There might be impor- 
tant nonperturbative features missed in every two-body confinement potential model. 
For example, three-gluon exchange between three quarks and three-body instanton in- 
teractions do not contribute within a colorless meson or baryon, but do contribute to 
a multiquark system P, 0. There are many other types of quark interactions which 
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result from multigluon exchanges which can not be included in a two body confinement 
potential. The QCD condensates within a nucleon and between two nucleons might be 
different. These arguments show that a direct extension of the two body confinement 
to the multiquark system has not been justified even though it might be a good ap- 
proximation in the single hadron case. The absence of experimentally observed color 
Van der Waals forces has been a problem for the two body confinement potential model 
for a long time. Multiquark systems provide more variations of low energy behavior 
of QCD and allow for further tests of the phenomenology built into quark models, 
especially the nature of confinement. 

The quark delocalization and color screening model (QDCSM) has been developed 
in which an unusual confinement parameterization has been introduced. Following 
the concept of electron delocalization in molecular orbitals, quark delocalization is 
introduced to enlarge the variational Hilbert space to include various deformed six- 
quark configurations. The usual pair of three-quark clusters and six-quark, "bag", 
configurations are the two extremes of this space |jlO| . 



With only one adjustable parameter, the color screening constant /i, the existing A^A^ 
(/J = 01, 10, 11, 00), A^A (/J = |0, il) and A^S (/J = |0, ^1, |0, |l) scattering data 
have all been fit qualitatively [p!0| , |ll], |12|, |13| within this model. A relativistic version of 
this model has recovered ^He, ^He, as approximate four- and three-nucleon systems 



starting from appropriate 12- and 9-quark configurations |T^, |15|, |T6 



Both the relativistic and the nonrelativistic version give almost the same dibaryon spec- 
trum, reproducing small deuteron binding and the zero energy A^A^ IJ = 10 resonance. 
An IJ^ = OS"*" state, called the d*, was found to be most interesting because of its large 
binding energy (~350 MeV) and small decay width (~1 MeV) fig, 0, |T8], |20[ . This 
state appears in and has been studied in many models. Kamae and Fujita obtained 
almost the same large binding [El] . Other models obtained smaller binding [B^, ESl E4 



The experimental situation regarding dibaryons remains unsettled. Precise np total 
cross section measurements provide a stringent limit on the A^A^ decay width of any 
dibaryon in the d* mass range . Similarly, the H particle has been studied in many 
models and searched for experimentally for more than 20 years without any indication 
of its existence p6|. There is an IJ^ = 00~ state, d', which seems to have experimental 
support ||2^ but which is hard to accommodate in a six-quark model space l2^, p9| . 



Newer experiments using simple systems have not confirmed the existence of a d' sig- 
nal |3^. A high- mass dibaryon state predicted by Lomon et al finds support from 
SATURNE pp scattering data |32|. 

This paper reports further study of the d* binding and decay width with the aim of 
providing more reliable information for experimental searches for this dibaryon. The 
paper is organized as follows: A brief review of the features of the QDCSM appears 
in Sect. 2. Our calculation method is presented in Sect. 3. Results and discussion are 
provided in Sect. 4. The last section gives a summary. 



2 



2 Model Hamiltonian and wave function 



The details of the QDCSM can be found in Ref.||10|, |T8[. Here only the model Hamil- 
tonian and wavefunction used in our calculations are repeated. 

The Hamiltonian for the 3-quark system is the same as the usual potential model [|1|, 
and for a 6-quark system, it is assumed to be 



Vt 



i=l i<j=l 

vr5(f) / 2^ ^\ ^ / 
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— ttcXi ■ Xj < l^g-Af^ 



if i,j occur in the same baryon orbit, 
if i,j occur in different baryon orbits. 



-fj and all other symbols have their usual meaning except the confinement 



where f = 

potential which will be explained below. 

It is well known that confinement is a nonperturbative QCD effect; in general one does 
not expect it to be described by a sum of two-body interactions. We model confinement 
as follows: In our approach, the fractional parentage expansion method is used and 
the matrix elements of the six-quark Hamiltonian are simplified to a four-body overlap 
and a two-body matrix element. We assume the following recipe to determine the two 
body matrix element of confinement: The interaction takes the normal, unscreened 
form (quadratic in r) when the interacting quarks always remain in the same baryon 
orbit (see Eq.(9) below), both before and after interaction; otherwise the interaction 
takes the screening form (second form shown in the last of Eqs.(|l])). Although this 
has not been demonstrated to be correct, it is more sophisticated than the usual, 
simple two-body confining interaction, and does include a physically reasonable model 
of non-local, nonperturbative effects. 

Even though we use a two-body interaction form to evaluate the matrix elements, our 
model is, nonetheless, not a potential model. It is an effective matrix element approach 
extended from bound states to scattering states. It does reduce to the usual two-body 
confinement interaction within a single hadron. It also has the usual meaning of a 
two-body interaction in the asymptotic regime although not in intermediate regions. 
The main physics introduced is the recognition that the confining interaction between 
two nucleons might be different from that within a nucleon. In particular, we represent 
the nonlocal, nonperturbative backflow of color (and pair creation, as seen in lattice 
studies [^) by screening of the confining potential. If this model assumption of the 
confinement is not a good representation of the physics involved, we would expect that 



3 



to appear in a disagreement between our calculational results and data. The quality 
of the fit to B-B scattering data referred to in the introduction shows that there is no 
obvious contradiction. Our predictions with respect to dibaryon states will provide a 
further test. 

We use the resonating group method (RGM) to carry out a dynamical calculation. 
Following the nomenclature of Ref.|2^, we write the conventional ansatz for the two- 
cluster wavef unction as 

1*6,) = [[^^B.^B.f''^'' ® Xl{R)] ' , (2) 

L 

where [a] = [222] gives the total color symmetry and all other symbols have their usual 
meanings. \1/ Bi is the 3-quark cluster wavefunction (after removal of the center of mass 
motion), 

/ 9 \3/4/ 9 n3/4 + I 

"^'-{i^) (i^) » ^" ")l..S:iBM. (3) 

where Xc{Bi) is the internal color wavefunction of the baryon and the Jacobi coordinates 
are defined as follows, 

Pi = r 1 - f2, P2 = f4 - fs, 

Ai = rg - -(fi + rs), Aa = fe - -{u + rg), (4) 
- 1 - 1 

RB^ = + + ^52 = 3(^4 + ^5+^6), 

R = R-Bi — Rb2^ Rc = 2^^Bi + Rb2)- 

From the variational principle, after variation with respect to the relative motion wave- 
function x(-R) = J2lXl{R), one obtains the RGM equation 

J H{R, R')x{R')dR' = E J N{R, R')x{R')dR', (5) 

where H{R, R'), N{R, R') are Hamiltonian and norm kernels, respectively. Their de- 
tailed expressions can be found in Ref. 



The energies, E, and the wavef unctions, x{R)i ^i-re obtained by solving the RGM equa- 
tion. In practice, it is not convenient to work with the RGM expressions. We introduce 
generator coordinates. Si, to expand the relative motion wavefunction, x{R)i 



L 

After the inclusion of the center of mass motion 



^c{Rc) = (^) e-i^^c, (7) 
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the ansatz, Eq.(0), can be rewritten as 



6q 



[m.sABi)m.sAB2)Y'Y'^{s 



/3=4 
n J 



where 4'a{Si),(f)[3{~Si) are the single particle orbital wavefunctions with different ref- 
erence centers, 



4>a{Si) 



1 



i^(r«-5,/2)2 



1 \ 3/4 



(9) 



With the reformulated ansatz, Eq.(||), the RGM equation Eq.(|) becomes an algebraic 
eigenvalue equation, 

Y.C,Mf = ET.C,.L'Nl^- (10) 
i,L j 

where iVjj , Hj^J^ are the Eq. (Q) wavefunction overlaps and Hamiltonian matrix ele- 
ments (without the summation over L'), respectively. By solving the generalized eigen 
problem, we obtain the energies of 6-quark system and corresponding wavefunctions. 

In the QDCSM, the single particle orbital wavefunctions are delocalized. To implement 
this here, we modify Eqs.(|^) as follows: 



'111 



MS^) ^ MSi,e) = [MS^) + eM-S^))/N{e), 
N{e) = \/l + e2 + 2ee-^'/4^'. 



It is straightforward to extend the method to multichannel coupling. In the multichan- 
nel case, the ansatz used is 



6q 



/flQ 3 6 



13=4 
J 



[XciB,)xc{B2)f^ 



(12) 



where k is the channel index. The eigen equation is similar to Eq.(|TO|), with an ad- 
ditional summation over k. For example, for IJ = 03, we have k = 1,2,3,4, corre- 
sponding to the channels AA = 3 L = 0, AA S" = 3 L = 2, NN S = 1 L = 2, and 
AA S =1 L = 2. 
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3 Calculation method 



To further simphfy the calculation of the matrix elements of the six-quark Hamiltonian, 
the physical bases, Eq.(H), are first expanded in terms of symmetry bases (group chain 
classification bases). Then the powerful fractional parentage expansion method is used 
to calculate the matrix elements between symmetry bases of the six-quark system. 
Finally, the matrix elements between physical bases are obtained by the transformation 
between the physical bases and symmetry bases. The details can be found in Ref.fS^, 

m- 

The partial width for d* decay into the A^A^ D-wave state is obtained using "Fermi's 
Golden Rule", 



' mt..mi. y'^'^i 



Mj,,Mj 



\ E ^,m,,^ml-Am% J \M\'dQ, (13) 



Mj.,Mi, 



where M denotes the nonrelativistic transition matrix element, and Mj. and Mj^ are 
the spin projections of the initial and final states. 

M={d*\Hj\[^r,,^M,Y'e'P-^), (14) 



where R is as above and p = h v^d* ~ 4m^ is the available relative momentum be- 



tween the nucleons as determined by the energy conserving (5-function in Eq. ([T3|) . The 
interaction Hamiltonian, Hi, here is Hj = Vr- 

By expanding the plane wave in terms of spherical harmonics, and taking into account 
angular momentum conservation, the transition matrix element can be put in the form, 

1 J \ 



M = -An{d*\Hj\Ci^^l^,,^ 



[^r,,^^,Y'Y\R){Y\p)r 



jLipR + h))- (15) 



where we have included the phase shift, 6l, induced by the interaction in the spherical 
Bessel function. Of course, the phase shift is only apparent at large separations. 

On the other hand, the relative motion wavefunction between two nucleons can be 
determined from either the RGM equation or the algebraic equation, Eq.(0). The 
relative motion wavefunction can be written as 

Xl{R) = V4^T.C^l[^) e-^^^'^'^h,{^^RS,)Y\R), (16) 



where iiiz) is the modified spherical Bessel function. If we can match the expansion 



E^-U^ e-^^^-''^^di,RS^ (IT) 
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to the function jiipR + ^l) for the separation in the calculated region, then the wave- 
function used in Eq.(|l5[) is just the wavefunction defined in Eq.(||) with an additional 
factor of '\/4ti(Y^{'p))*, which contributes a factor An to the decay width after the 
integration over the angular space, Q. 

We have used a type of box normalization method to calculate the transition matrix 
element. Because the wavefunction of the initial bound state d* is compact, its ampli- 
tude at large separation (i? > i?o ~ 3.3 fm) can be safely neglected. By contrast, the 
final NN D-wave state with the same energy as that of the d* is a scattering state, 
so it extends over the entire space. However, the transition matrix element depends 
only on the overlap between the d* and A^A^ D-wave states, so only that part of the 
wavefunction of the NN D-wave with separation less than i?o fm can contribute. To 
determine the wavefunction of the final state NN D-wave in this finite range, the RGM 
equation is used so that the final state interaction is taken into account automatically. 
By adjusting the box size, the correct internal relative momentum A^A^ D-wave state 
can be obtained to meet the requirement of energy conservation in the decay. Finally, 
a normalization correction due to the finite range of the wavefunction is needed, which 
is obtained by matching the expression Eq.(|l^ with jiipR + ^l)- In this way, we are 
able to address the bound state problem and the scattering state problem on the same 
footing. For example, the phase shifts of the scattering state can be obtained from the 
comparison of the calculated wavefunction with the free Bessel function, 

4 Results and discussion 

As a test of the model, we first attempted a calculation of the deuteron. The parameters 
we used in the calculation, which are fixed by baryon properties and A^A^ scattering, 
are 

m„ = mrf = 313 MeV, h = 0.603 fm, a = 25.13 MeV/fm^, a, = 1.54, ^ = 1.0 fm'^ 

These parameter values are very similar to those usually appearing in constituent 
quark models [Q]. The single channel and channel coupling calculations do show an 
attraction between nucleons, but it is insufficient to form a bound state. The mixing 
of the D-wave into the deuteron state is also too small (less than 1%) to affect the 
result. The deuteron can be forced to form in the QDCSM by increasing the color 
screening parameter. However the small radius of the resulting object and negligible 
D-wave mixing demonstrate that it is not physically correct to do this. 

Together with the absence of a long range tail for the A^A^ interaction, this serves 
to emphasize the need for the long range part of one pion exchange (OPE), which is 
absent from our model. To reproduce the deuteron, which is an extremely extended 
object, it will be necessary to extend the QDCSM to include OPE. To avoid double 
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counting, a short- distance cutoff of the pion-quark couphng is also required. Prehmi- 
nary calculations show that the deuteron can be well reproduced in this manner. We 
will address these issues in a future publication p6 |. 

The case of the d* is quite different, however. The small size and the large delocalization 
of quarks shows that it is a highly compact object. The effect of OPE should be 
considerably reduced compared to the case of the deuteron. Using the same parameters, 
a dynamical calculation of the d* was carried out. These results appear in Table I. We 
find that the mass of d* is somewhat larger than that obtained by our earlier adiabatic 
calculation |]l8l, mainly because the relative motion energy found is larger than that 
of the zero-point oscillation energy used in the adiabatic calculation. Although this d* 
mass exceeds the NNinr threshold, it should remain a narrow resonance due to the 



extremely small phase space available [B7 



sc= 



Table I. Mass and radius of the d* in the QDCSM. 
single channel; cc2 = AA S-wave + NN D-wave; cc4 = full four-channel 

coupling. 





sc 


cc2 


cc4 


adiabatic [|18| 




mass (MeV) 


2186 


2180 


2176 


2134 


V< > (fm) 


1.2 


1.2 


1.3 


1.4 



We have also calculated the partial decay width of d* into the NN D-wave, using the 
method introduced in last section. The result is dominated by the stretched (mj = ±3) 
states which contribute P^^^' = 1.95 MeV each to the sum over spin projections. We 
find P = 4.32 MeV for the total width. This result is similar to that obtained from 
meson exchange pO[. It should be expected that the total width will be increased by in- 



cluding OPE as discussed above for the deuteron. Although it is generally believed that 
the tensor interaction due to OPE is much stronger than that from one gluon exchange, 
we expect this increase will be limited due to the combination of the compactness of 
the d* and the required short- distance cutoff of the pion-quark coupling. 



5 Summary 

In the framework of QDCSM, dynamical calculations of the d* dibaryon and the partial 
decay width of the d* into an NN D-wave state were carried out. The results include a 
slight increase in the mass of the d* over the adiabatic result, rising above the NNtttt 
threshold. However, the small phase space still suggests that the resonance will be 
narrow and, indeed, the calculated d* — > A^A^ D-wave decay width is on the order of a 
few MeV. 

Our calculation of the deuteron shows the predictive power of the QDCSM on the one 
hand, but on the other, that OPE must be included in the calculation of extended 
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objects (with small delocalization). However, we do not expect that the addition of 
OPE will have a large effect on the calculation of d*, because of its compactness. 
Addition of OPE to the QDCSM is underway and preliminary results show that both 
the D-wave component and the radius of the deuteron are indeed improved. Details 



will be presented in a forthcoming paper |36 



The calculated d* NN D-wave decay width is larger than allowed by the Lisowski [25 
data. However for a state as high as 2.2 GeV, sea quark excitation or the NNir com- 
ponent may well be important. Inclusion of this component may affect the mass and 
decay width strongly. In particular, the d* — > NNir decay width may be significantly 
increased [0. We believe that it is still too early to conclude that the d* has been 
ruled out by precise np total cross section measurements and that further calculations 
including these complications are needed. 
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